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Abstract. The inclusive production cross sections of the charmed mesons D°,D+,D+ and 
D* + have been measured in interactions of 920 GeV protons on C, Ti, and W targets with 
the HERA-B detector at the HERA storage ring. Differential cross sections as a function 
of transverse momentum and Feynman's x variable are given for the central rapidity region 
and for transverse momenta up to pt = 3.5 GeV/c. The atomic mass number dependence 
and the leading to non-leading particle production asymmetries are presented as well. 



1 Introduction 

The cross sections for charm and beauty hadro-production are of considerable theoretical interest 
[TJEJE] ■ Perturbative QCD is expected to work well for the large mass top quark production and 
less well for the lower mass b and c quarks [H0E]. At present, several published results of 
measurements of charm production in proton- nucleus collisions [4,5,6,7,8,9] are available. They 
are mainly restricted to beam energies between 200 GeV and 800 GeV and mostly have low 
statistics. Only one of the experiments also provides a measurement of the dependence of the 
cross section on the atomic mass number. More data would help in determining the strong 
interaction parameters as well as in guiding the calculation of non-perturbative effects. Another 
motivation comes from the prediction that a prominent manifestation of the quark gluon plasma 
at the LHC is a larger ratio of charmonium to open charm cross section compared to, e.g., 
production in pA collisions at lower energies [10J. The present work provides a new data point 
at 920 GeV proton beam energy. 

Collisions of the 920 GeV HERA proton beam in C, Ti and W fixed targets have been 
measured with the HERA-B spectrometer. In previous papers we have reported on the bb 
T [12], and charmonium [131114] production cross sections, while the present work deals with 
production of open charm in the inclusive reactions pA — > DX. Here D represents a D°,D + , Df 
or D* + detected through the decay channels: D° — ► K~ir + , D + — > K~tt + it + , Df — > 07r + — ► 
K~ K + 7r + , and D* + — ► D°7t + — > K~ir + ir + . Throughout this paper, charge-conjugated modes 
are included unless noted otherwise. 

The paper is organized as follows. We first briefly describe the apparatus, the data sample 
and the method of analysis. We then present the results and finally make a comparison with 
other measurements and theoretical expectations. 



2 The detector 

HERA-B was a fixed target spectrometer (see Fig. [T]) using the 920 GeV proton beam of the 
HERA e-p collider. Interactions occurred on one or more wires (depending on run configuration) 
which were organized into two groups of four target wires each; the groups were separated by 
4 cm along the beam, and the transverse wire dimension was 50 /im-500 fim [15]. The wires were 
positioned in the beam halo, and their distance to the beam core was automatically adjusted 
to maintain a constant interaction rate. Details of the various subdetectors have been published 
[T6 l[T7][18| ll9 y20ll2T] , so only a brief overview of the apparatus is given here. 

Tracks originating from proton interactions and decay vertices were measured with a vertex 
detector system (VDS) [16] . Sixty-four silicon strip detectors (50 x 70 mm 2 , pitch of ~50 //m) 
with double-sided readout were arranged in eight stations between 7 cm and 200 cm downstream 
of the targets. The detectors were in Roman pots [22] under vacuum and their inner edges 
were adjusted to be in the range 10-15 mm from the beam center. With this system, a vertex 
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resolution of a z ~ 500 along the beam direction and a x<y ~ 50 fxm in the transverse plane 
was achieved. 

Particle momenta were measured with a tracking system and a dipole magnet of 2.13 T-m 
field integral. The first tracking chamber was upstream of the magnet, and the remaining six 
chambers were downstream of the magnet between 7 m and 13 m from the interaction region. 
Due to a large variation of particle flux density, the tracking system was divided into a fine- 
grained inner tracker (ITR), using microstrip gas chambers with gas electron multipliers and 
~300 fim pitch [IT], and a coarse-grained outer tracker (OTR), using honeycomb drift cells with 
5 mm and 10 mm cell diameters [18] . The obtained momentum resolution can be parameterised 
as (Tp/p = (1.61 + 0.0051 • p [GeV/c]) % [18] , where p is the particle momentum. 

Particle identification was achieved with three subdetectors: a ring imaging Cherenkov counter 
(RICH), an electromagnetic calorimeter (ECAL), and muon chambers (MUON). The RICH de- 
tector [19] is a large vessel containing about 100 m 3 of C4F10 gas at STP, which provided about 
2 m of radiation path. The Cherenkov photons were focused and reflected by two sets of spherical 
and planar mirrors onto an upper and a lower photon detector, located well outside of the main 
particle flux. Each photon detector consisted of about 1100 multi-anode photomultiplier tubes. 
The identification efficiency for pions was about 90% in the momentum range from slightly above 
the pion threshold (2.4 GeV/c) up to 70 GeV/c, with the kaon misidentification probability al- 
ways below 10%. The efficiency for kaons was above 85% for momenta between 15 GeV/c and 
45 GeV/c, with a pion misidentification probability of ~1%. For particles below the Cherenkov 
threshold of ~10 GeV/c, the misidentification of pions as kaons was kept below 10%. 

The ECAL [20] was a sampling calorimeter of the "shashlik" type, with scintillator plates 
sandwiched between tungsten (inner region) or lead (outer region) absorbers. The calorimeter 
was read out by optical fibers and photomultiplier tubes with the readout granularity of inner and 
outer regions adapted to different particle rates in order to maintain acceptable occupancies. The 
MUON system [21], situated in the most downstream region, consisted of four large detector 
stations separated by concrete and iron absorbers. Each muon detector plane had gas pixel 
chambers in the inner region and gas proportional tubes in the outer region. 

For the present measurement, the analysis of data was based on the vertex detector, the OTR 
tracking system and the RICH counter. 
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Table 1. Summary of the data statistics and the integrated luminosities of the present study. 



Target 


A 


events [xlO 6 ] 


C [Mb" 1 ] 


C 


12.01 


89.3 


375 


Ti 


47.88 


24.7 


31 


W 


183.84 


67.6 


36 



3 Data analysis 

The analysis was performed on data sets with a single target wire made either of carbon, titanium 
or tungsten. Only runs with stable conditions and a minimum bias trigger were considered; the 
resulting sample consisted of 182 million interactions (Table [T]). The trigger required at least 20 
hits in the RICH detector (compared to an average of 33 for a full ring from a (3 = 1 particle 
|19j ) and had an efficiency etriggcr ~ 95% for inelastic interactions. The integrated luminosity 
L was determined |23| from the number of measured inelastic interactions iV; ne i and the total 
inelastic cross section <7i ne i, using the expression C = iV ine i/(et r igger<7inei)- The data were recorded 
at a moderate interaction rate of about 1.5 MHz which corresponds to 0.17 interactions per filled 
bunch crossing. Therefore only about 10% of triggered events contain more than one interaction. 
The data acquisition rate was about 1 kHz, and the bulk of the data was recorded within a two- 
week running period. 

At the HERA-B energy, the charm production cross section is more than two orders of magni- 
tude smaller than the inelastic cross section. Taking into account the relatively small branching 
ratios for the D meson decay modes into two or three charged particles, one expects sizeable 
backgrounds. Particle identification alone is not sufficient to extract signal events. However, the 
large boost of the center-of-mass system of HERA-B (7 = 22), causing D mesons to decay sev- 
eral millimeters from the target, combined with good vertex resolution (~0.5 mm longitudinally) 
allows us to distinguish D meson decay products from particles originating at the primary in- 
teraction point. The data selection thus requires a detached secondary vertex formed by tracks 
not coming from the primary interaction point as well as the identification of kaons and pions. 
The selection criteria are summarized in Table [2] and discussed in some detail below. 



3.1 Data selection 

In addition to a detached vertex for the ground state D mesons, at least one reconstructed 
primary vertex was required in each selected event. Primary vertices were determined from all 
track segments reconstructed in the VDS. Since the proton interaction point must be inside the 
target wire, the primary vertex coordinate transverse to the wire direction was replaced with 
the known target wire position. 

For the tracks corresponding to the decay products of D mesons, mild requirements were 
applied on the number of hits in the vertex detector and in the main tracking system, as well as 
on the track fit quality. Selection criteria used to identify the final state kaons and pions were 
based on information from the RICH counter. While a rather strict cut on the kaon likelihoocQ 
was required for the kaon candidates (Lk > 0.5 for the D° and D + selection, and Lk > 0.33 
for Df candidates), for pions only a mild cut was applied on the sum of RICH likelihoods for 
electrons, muons and pions, L e + + L n > 0.05. No particle identification requirement was 

Likelihoods for the electron, muon, pion, kaon, pro- 
ton and background hypotheses are normalized such 
that their sum is equal to one: L e + L M + + Lk + 

Lp + -/vbkg = I- 
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imposed for pions from D* + — > D tt + decays, which tend to have low momenta and are thus 
denoted as vr s i ow . 

The tracks were combined to form D°, D + and Df candidates. Candidates with an invariant 
mass in the interval of ±500 MeV/c 2 around the D meson nominal mass were retained for 
further analysis. For the Df candidates, the invariant mass of the K + K~ pairs was required to 
be in the interval of ±10 MeV/c 2 around the <p nominal mass; the absolute value of the cosine 
of the angle 9^ between the K + and tt + in the rest frame of the was restricted to the values 
above 0.5, exploiting the vector nature of the intermediate state (p. The D* + candidates were 
reconstructed from D° candidates with invariant mass within ±75 MeV/c 2 (3.5a in resolution) of 
the D° mass and slow pion candidates, after a vertex fit to the D°; an additional cut was applied 
on the product of transverse momenta of the D° daughter tracks and of the D° momentum, 
p(D°) PT (K)p T (iT). 

In addition to the criteria described above, the analysis was restricted to the region of phase 
space with high acceptance, —0.15 < xp < 0.05, where xp is the Feynman x variable. The 
daughter tracks of a D meson candidate were fitted to a common vertex. Only combinations 
with a vertex probability greater than 0.1% and with a secondary vertex displaced by more than 

4 standard deviations downstream of the wire were accepted. 

The D meson candidate was then associated with the primary vertex. In the case of events 
with multiple reconstructed primary vertices, the vertex with the smallest impact parameter 
significance (i.e., the measured value divided by its estimated error given by the covariance 
matrix) with respect to the track of the D meson candidate was chosen. To avoid a possible bias 
in the primary vertex position due to tracks from the D meson candidate, the primary vertex 
was re-fitted without the D daughter tracks. 

The final set of criteria was based on the primary and secondary vertices. The main source 
of background, which is due to combinations of particles emerging from the primary interaction 
point, was reduced by the following requirements: (1) the secondary vertex should be detached, 
(2) the tracks forming the secondary vertex should not come from the primary interaction point 
and (3) the D meson candidate should originate from the primary interaction point. To fulfill 
these criteria, cuts were applied on the following variables: 

— d(D) the significance of the distance between the secondary vertex and the associated primary 
vertex, 

— b(Tr),b(K) the significance of the impact parameter of a pion or kaon with respect to the 
primary vertex (in case more than one primary vertex was found, that closest to the particle 
is chosen), 

— 6(D) the significance of the impact parameter of a D meson candidate with respect to the 
associated primary vertex. 

We found that, for the three-body decays of D + and Df, a more effective cut than a selection 
based on single impact parameter significances of daughter tracks is a cut on their product 
b(K)b(-Ki)b(ir2) and b(Ki)b(K2)b(-7r) , respectively. The background level is further reduced with 
criteria of the form ^6(^)6(^1)6(^2) > 4(t - t ), for D+, and v /6(^)6(tt) > 0.75(i - to), for D+. 
Here the proper lifetime t is in units of the D mean lifetime, and the offset to is determined in 
an optimization. 

For each decay mode the optimal cuts were determined by maximizing the signal significance 
S/y/S + B in a ±3<7 window centered at the D meson nominal mass (signal window). The signal 

5 was taken from Monte Carlo simulation and was scaled to the luminosity of real data by using 
an estimation for the production cross sections from fits to the published D meson cross sections 
PP. For the the cross section was assumed to amount to 20% of the sum of the D° and D + 
cross sections [2]. 
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In the case of ground state D mesons, the number of background events B was estimated 
from the data sidebands. For the D+, the mass region of ±50 MeV/c 2 around the nominal mass 
of the D + , where a contribution of the decay D + — ► (J)tt + is expected, was excluded from the 
lower sideband. To reduce the sensitivity to statistical fluctuations, the sidebands were chosen 
to be larger than the signal window. In the case of the D* + , the wrong sign combinations from 
real data were used to estimate the background in the signal window. As is usually done for this 
decay, the signal was reconstructed via the mass difference q = m(K,TT,ir s \ ow ) — m(K,Tr) — m n 
rather than using the invariant mass of the K, tt, ir s \ ow combinations. 

Table 2. Selection criteria; d and b denote decay distance and impact parameter significances, respectively, and 
Lh are likelihoods for a hypothesis h. 





D+ -» K--K+-K+ 


chosen a priori 


L K {K) > 0.5 

L e (w) + L /J ,(n) +L n (n) > 0.05 


L K {K) > 0.5 

L e (%) + L m (tt) + L^(tt) > 0.05 


optimized using background data and signal MC 


d{D°) > 6.1 
b(D°) < 2.4 
b{K) > 3.4 
&(tt) > 3.7 


b(D+) < 2.6 
6(^)6(^)6(^2) > 106 
^6(if)6(7ri)6(7r 2 ) > 4(t - to), to = 2.48 








chosen a priori 


L K {K) > 0.33 

L e (ir) +Mtt) +L^(tv) > 0.05 
\m{K+K-) - m{cj>)\ < 10 MeV/c 2 
cos(^)| > 0.5 


L K (J?) > 0.5 

L e (Ti) + L m (tt) + L„(n) > 0.05 
Im(A'Tr) - m(D°)\ < 75 MeV/c 2 
d(D°) > 4 


optimized using background data and signal MC 


dipt) > 5.3 
b(Dt) < 2.11 
b(K-)b(K+)b(iv) > 28.7 
^/b(<t>)b(ir) > 0.75(t - t ), to = 1.0 


6(D°) < 2.4 
b(K) > 2.1 
6(tt) > 1.7 

p(D°)pt(K) P t(ty) > 17.7(GeV/c) 3 



After applying the selection criteria, which are summarized in Table [21 the remaining data 
were scanned for events with more than one D meson candidate (in 0.5%, 8%, 0% and 20% 
of events for the D°, D + , Df and D* + candidates, respectively). In case of the D°, D + and 
candidates, the combinations with the largest decay distance significance d{D) were kept. 
For the D* + analysis, first the candidates whose intermediate D° had the largest decay distance 
significance were selected. If multiple candidates remained (i.e., due to multiple vr s i ow candidates), 
the D* + candidate with the highest vertex probability was kept. 

3.2 Signal yields 

The invariant mass distributions for D meson candidates after applying the selection criteria 
discussed above are shown in Figs. [20 

The signal yields are extracted from the histograms by a maximum likelihood fit assuming 
Poisson statistics in individual bins. A Gaussian function is used for the signal, while the back- 
ground description depends on the type of the D meson. For the D + and Dj, the background 
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is fitted by an exponential function. The background for D candidates is more complex and 
consists of a combinatorial part, fitted by an exponential, and a contribution from partially 
reconstructed charm decays. This latter background is visible in the mass range below the D° 
peak. The shape of this background is taken from Monte Carlo simulation of cc events to which 
the same selection criteria are applied as for the data. 

In the Df invariant mass distribution the Cabibbo suppressed decay of D + — > <pir + is also seen 
(the peak to the left of the Df peak). This peak is included in the fit function as an additional 
Gaussian of the same width as the signal Gaussian, with its normalization as an additional 
free parameter, and its mean fixed to that extracted from the D + — > K~tt + tt + invariant mass 
distribution. The background for the D* + candidates is parameterised as a(g 1//2 + bq 3 ^ 2 ), with 
a and b as free parameters. 

The fitted peak positions are within one standard deviation of the corresponding world aver- 
age values [24] . with the exception of D+, which deviates by less than two standard deviations. 
The widths of the signal peaks are about 30% larger than the corresponding Monte Carlo values. 

The numbers of reconstructed D mesons are summarized in TableEl In total, 175 D°, 131 D + , 
11 Df and 61 D* + decays are found. The yields for various subsamples (particle, anti-particle, 
individual target material) are obtained by fitting with the mean and the width (r.m.s.) of the 
signal Gaussian function fixed to the values obtained from the fit to the full sample. 




1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 

m(Kn) [GeV/c 2 ] 




1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 



m(Klz%) [GeV/c 2 ] 

Fig. 2. Invariant mass distributions for K~ir + (a) and K~tt + tt + (b) combinations. The curves show results of 
maximum likelihood fits to the data. 



The D meson proper time distributions provide a check for signal consistency. The acceptance- 
corrected distributions are found to be consistent with the expected exponential decay. To de- 
termine the lifetimes, a simultaneous likelihood fit of events in the signal window and in the 
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D 



m(KKn) [GeV/c] 

Fig. 3. Invariant mass distributions for (f>n + — > (K~ K + )tt + combinations. Besides the Df peak at 1.96 GeV/c 2 , 
a D + peak at 1.87 GeV/c 2 is also visible. This peak corresponds to the Cabibbo suppressed decay D + — * (f>ir + . 
The curve shows the result of a maximum likelihood fit to the data. 




2.5 5 7.5 10 12.5 15 17.5 20 22.5 25 



q [MeV/c 2 ] 

Fig. 4. Invariant mass difference q — m(K, ir, 7r s i ow ) — m(K, ir) — m w for K~ir + ir + combinations. The curve shows 
the result of a maximum likelihood fit to the data. 

Table 3. Number of reconstructed D mesons. 



sample D° D + DJ D 



total 


174.8±16.8 


130.5±14.7 


11.4±4.0 


61.3±13.0 


particle 


75.9±10.9 


54.5±9.3 


4.9±2.6 


21.0±6.6 


anti-particle 


99.0±11.9 


75.8±10.5 


6.7±2.8 


40.6±8.3 


C 


66.1±9.6 


43.1±7.7 


4.2±2.2 


26.6±6.4 


W 


92.3±11.7 


72.4±10.6 


6.7±3.0 


24.8±7.5 


Ti 


17.4±5.7 


14.9±5.0 


0.4±1.0 


9.6±4.0 



sidebands is used. The results, cr = (302±33) fim for D + , (120±13) for D°, and (165±52) fim 
for Df, are in good agreement with the world average values [2~4*] . 

As a further consistency check of the and D* + signals, the intermediate states eft and 
are checked. They should be visible in the corresponding invariant mass distributions, when 
the signal region in the initial state invariant mass distribution is selected, and the fit to the 
intermediate state invariant mass distribution should give an event yield consistent with the 
number of events in the initial state peak. The yields extracted in this way are in good agreement 
with the values given in Table [3l For the Df decays, the distribution of the cosine of the angle 
between K + and tt in the rest frame of the (j) is found to be consistent with the expectation 
given by the vector nature of the intermediate state 4>. The expected D + — > cj)7r + signal yield 
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can be estimated from the number of reconstructed decays from our measured D + cross section 
in the D + —* K~ir + Tr + decay channel. The fitted number of events, 9.8±3.8, is in reasonable 
agreement (1.4<t higher) with the estimated number, (4.2±1.2). 



4 Efficiency determination 

A Monte Carlo simulation is used to determine the signal reconstruction efficiencies. The Monte 
Carlo samples for pA — ► DX are generated in two steps. First, a cc pair is generated with 
Pythia 5.7 [25] such that a particular D meson is always produced. The generated events are 
re-weighted to make the resulting cross sections conform to the parameterisations 



da 
dp\ 



1 + 



l)PT 



2 r(p - 1) <pt> 



and 



da 



Aexp (- 



dx F I A'(l 



M) 



\Xf 

\xf 



< x h , 
> x h , 



(2) 



with o n 



and In £ 



'g ~ \J n awj> - 111 W ~ n ^ 2(i-x b ) + m (l — x b)] The average transverse momentum 

of (px) = 1-04 ± 0.04 GeV/c and the value of the exponent /? = 7.0 ± 4.3 are taken from the 
present analysis (Sec. I5.2p . The value n = 7.7 ± 1.4 is taken from the average of E743 and E653 
[BJ[7], and x\, is assumed to be 0.062 ± 0.013 as measured by the E791 experiment [26J. The 
influence of the parameter uncertainties and other possible parameterisations of Eqs. [TJ and [2] 
are taken into account as systematic errors. After the generation of the D mesons, the remaining 
energy is input to the Fritiof 7.02 [27] program package which generates the underlying event 
taking into account further interactions inside the nucleus. 

The detector response is simulated with the Geant 3.21 package J2SJ . Realistic detector ef- 
ficiencies, readout noise and dead channels are taken into account. The simulated events are 
processed by the same reconstruction codes used for the data. 



5 Results 

5.1 Total cross sections 

The visible cross section per nucleus, i.e., the cross section measured in our visible range of 
—0.15 < xf < 0.05, is given by 

A(7pA > 1 = Br-*- (3) 

where Ni is the number of reconstructed D mesons for a particular target i, e« and Ci are the 
corresponding efficiency and integrated luminosity, and Br is the world average branching ratio 
for a specific decay channel [23]. The cross section for D meson production on a nuclear target 
of atomic mass number A is parameterised as 



cr p A = c P N • A a . 



(4) 
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Table 4. Summary of systematic uncertainties of visible cross sections 



Source D° D+ D+ D*+ 

Event fitting 3.4% 2.6% 6.0% 9.7% 

Branching fractions 1.8% 3.6% 13.0% 1.9% 

Luminosity 3.7% 3.7% 3.7% 3.7% 

Reconstruction efficiency 8.4% 10.3% 12.9% 9.7% 

Total 10% 12% 20% 14% 



To combine data recorded with different target materials, the production cross sections per 
nucleon Z\<t p n are extracted in the following way. From Eq. [3] and Eq. HJ the D meson yield of 
the target i is derived, 

Ni = Br • e; • d ■ Z\a pN ■ Af. (5) 
By summing Eq. E] over all targets and solving it for the production cross section we get 

N 

where N = J2i N is the measured D meson yield of the total data sample. The sum in the 
denominator of Eq. [6] can be rewritten by introducing the average efficiency e, defined by the 
weighted sum 



Then the expression E] reads 



^Wi€i , Wi= 1 . (7) 



N 

^ N = Br-e-E^f (8) 



Since there is no experimental indication for nuclear effects in open charm production, a linear A 
dependence of the cross sections is assumed, and a is set to one in Eqs.[T|and[8l This assumption 
is also in agreement with our measurements as discussed in Section [5.51 

The total systematic uncertainties are, according to Eq. El composed of contributions from 
uncertainty in the signal yields associated with the fitting procedure (listed as 'event fitting'), 
branching fractions, integrated luminosity and reconstruction efficiency. The uncertainty of the 
reconstruction efficiency can be further divided into contributions from Monte Carlo statistics, 
track reconstruction efficiency, particle identification efficiency, selection criteria and the contri- 
bution from the re- weighting of kinematical distributions based on and xp. The individual 
contributions are summarized in Tables HHSJ The uncertainty due to selection criteria is de- 
termined in two ways: by varying the cut values in the selection criteria discussed above, and 
by performing a second analysis with an independent set of selection criteria. The change in 
the resulting cross section is taken as the corresponding systematic uncertainty. The systematic 
errors of track reconstruction and RICH particle identification (1.5% and 2.0% per track corre- 
spondingly) are estimated using decays of Kg and <fi as sources of pions and kaons. Note that 
no systematic error is assigned to the assumption a = 1 in order to be compatible with previous 
experiments. Using Eq. Eland data in Tables Q] and El Aa p ^ for a ^ 1 can be extracted. 

The resulting cross sections in the visible range, Aa P N and Aa v A, are summarized in Tables El- 
[HI In order to extrapolate the measurements to the full phase space, 

CTpN = — ( , (9) 
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Table 5. Summary of systematic uncertainties of reconstruction efficiency. 



Source 


D° 


D+ 


Df 


D*+ 


Monte Carlo statistics 


1.2% 


1.3% 


3.7% 


1.1% 


Track reconstruction 


3.0% 


4.5% 


4.5% 


4.5% 


Particle identification 


4.0% 


6.0% 


6.0% 


4.0% 


Selection criteria 


6.0% 


6.0% 


6.0% 


6.0% 


Re- weighting 


2.9% 


3.6% 


7.7% 


4.6% 


Total 


8.4% 


10.3% 


12.9% 


9.7% 



the fraction / v i s of D mesons in the visible range, defined by —0.15 < xp < 0.05, is determined 
in the following way. With the value of the exponent n = 7.7 ± 1.4 as measured by the E743 and 
E653 experiments [ElE], / v i s is calculated by using Eq. El The values corresponding to the choice 
x h = 0.062 [26] and x h = are 0.542±0.048 and 0.558±0.051, respectively. The difference is 
small compared to the uncertainty due to the error of the parameter n. For the extrapolation the 
average of both numbers, / v i s = 0.55 ± 0.05, is used. The resulting cross sections extrapolated 
to the full phase space are listed in Table EJ 

Figure [5] shows the comparisons of these measured cross sections with other experimental 
studies. To compare the data points at different center-of-mass energies, the overall normalization 
of theoretical predictions given in p] is determined separately for each of the D mesons by fits to 
the data which also include the results of the present measurements p. Our results are consistent 
with previous studies and represent an improvement at high energies. Note that our a{D + ) is 
somewhat lower than expected, leading to a lower ratio a(D + ) /a(D°) as will be discussed below. 

The measured sum of the D meson cross sections per nucleon a(D°) + a{D + ) + a(D+) = 
(87.4±8.2± 12.6) fib is used to determine the charm production cross section. From the fractions 
of charm hadrons produced in the hadronization of c quarks as measured in e + e~ collisions, the 

2 The original measured values are rescaled in [T] by 
taking into account the updated D branching fractions. 
The result of a(D°) by E789 is excluded from the fit. 

Table 6. Cross sections in the visible range (—0.15 < xf < 0.05) and extrapolated to the full phase space. 
The first error is statistical and the second systematic. In the second column, the systematic error due to the 
extrapolation uncertainty is quoted separately. 





ZW P N[/ib] (visible range) 


cr p N[/ib] (full if range) 


D° 


26.8±2.6±2.7 


48.7±4.7±4.9±4.4 


D+ 


11.1±1.2±1.3 


20.2±2.2±2.4±1.8 


Df 


10.2±3.5±2.0 


18.5±6.4±3.7±1.7 


D*+ 


11.9±2.6±1.7 


21.6±4.7±3.0±2.0 



Table 7. Cross sections for particle and anti-particle production in the visible range (—0.15 < ip < 0.05). The 
first error is statistical and the second systematic. 







particles 


anti-particles 


D° 


12.0±1.7±1.2 


14.8±1.7±1.5 


D+ 


4.8±0.8±0.6 


6.3±0.9±0.8 


Df 


4.1±2.3±0.8 


6.3±2.6±1.3 


D*+ 


4.5±1.4±0.6 


7.2±1.5±1.0 
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Table 8. Visible cross sections per nucleus. Numbers in parentheses in the last row are for the subsamples of D* + 
with D° daughters not common to the D° samples of the first row. The first error is statistical and the second 
systematic. 



Aa p A [mb] 



data sample 


C 


Ti 


W 


D° 


0.36± 0.05± 0.04 


1.01± 0.33± 0.10 


4.91± 0.62± 0.54 


D+ 


0.12± 0.02± 0.02 


0.48± 0.16± 0.06 


2.17± 0.32± 0.28 


Df 


0.12± 0.07± 0.02 


0.15± 0.35± 0.03 


2.18± 0.99± 0.46 


D*+ 


0.17± 0.05± 0.03 


0.71± 0.30± 0.10 


1.71± 0.51± 0.26 




(0.21± 0.06± 0.03) 


(0.95± 0.50± 0.14) 


(2.11± 0.84± 0.34) 
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Fig. 5. Comparison of the present results with previous measurements 4 5,6,7,8,9 . In figure a) the result of 
E789 is excluded from the fit. The different curves correspond to different assumptions on parton distribution 
functions [T]. 



production of D , D + and D s mesons is found to account for fjj = 0.891 db 0.041 of the charm 
cross section |29| . where the uncertainty in fjj is obtained from errors in individual fractions 
neglecting possible correlations. Assuming the same fraction, the present study derives the charm 
cross section per nucleon, a(cc) = (a(D°) +a(D + ) +a(Df))/(2f£>), where the factor 2 accounts 
for the charge-conjugated states which are included in the D production cross sections. The 
resulting charm cross section per nucleon is thus a(cc) = (49.1 ± 4.6 ± 7.4) ^b. Note that due 
to correlations, the systematic error in the sum of cross sections is larger than the value which 
one would get by adding in quadrature the individual contributions. 
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5.2 Differential cross sections 

The differential cross sections do jdp\ and da/dxp are determined from production yields in 
bins of p\ and xp by using Eq. [HJ The yield in each individual bin is determined by subtracting 
from the number of events in the D meson signal window the number of background events as 
estimated from the sidebands. The resulting differential cross sections for the production of D° 
or D + mesons are shown in Fig. [6J 

The parameters of the measured differential cross sections are determined in the following 
way. Because of the low statistics in individual bins of p\ and xp, we do not directly fit the 
background-subtracted distributions but instead do simultaneous binned likelihood fits of the p\ 
and xp distributions of events in the mass signal and sideband windows. This correctly accounts 
for the Poisson errors. 

For the transverse momentum distribution, several parameterisations can been found in the 
literature [2ll26tl3U]. In the present analysis, the parameterisation given in Eq. [T] is used, since 
our previous studies of J ftp, K* and (j) production indicate that it fits the data well over a large 
range of p\. The values of parameters (pr) and (3 are extracted from the fit. The distribution of 
background events Nt 0& -{p\) are assumed to be the same in shape and normalization for the signal 
window and sidebands. Several parameterisations of the distributions of events in sidebands fit 
well. We use the parameterization with two free parameters which gives the smallest x 2 / n df: 
Ce Bp Tior D° and D+. The same parameterisation multiplied by the efficiency in p\, 
is used for D* + . 

The D°, D + and D* + data samples are fit simultaneously, where in the D* + case only the 
subsample with D° daughters not common to the D° sample is included. The resulting fit 
parameters are (px) = (1-04 ±0.04) GeV/c and [3 = 7.0 ±4.3, with x 2 /ndf=0.86. The measured 
(pt) is significantly larger (by 3.3<r) than the value extracted from the Pythia Monte Carlo 
samples (0.90 GeV/c), while (3 is, within one standard deviation, equal to the value extracted 
from the simulation {(3 = 4.80). Fitting with the fixed value of (3 = 6, i.e., the value used for our 
study of J/ip production [13] yields the same (pt) value, and the same x 2 / n df. 

The xp distribution is usually parameterised with a power-law function: 

>hT oc(l-|x F |) n . (10) 



dxp 

This function differs in the central region from the predictions of the next-to-leading order QCD 
calculations [HJ. The measurements made by E791 [55J in 500 GeV irA collisions using a high 
statistics sample of 80k reconstructed D° mesons also show a similar discrepancy. They obtain 
an improved fit with the function given in Eq. 0, which uses a power-law function in the tail 
region and a Gaussian in the central regiorH. 

In our analysis of the measured differential cross section da/dxp (Fig. [5Jb)), the boundary 
parameter Xb was fixed to the value Xb = 0.062 as measured by E791 [55] . because our range 
of —0.15 < xp < 0.05 is too small to determine this parameter. The exponent n, is extracted 
by simultaneously fitting the D° and D + data samples with a likelihood fit to the events in the 
signal windows and sidebands. While the fitted value of the exponent n = 7.5 ± 3.2 agrees with 
the results of E653 [7] and E743 [6j , the statistical error is larger in our study. 

5.3 Ratios of cross sections 

The measured cross section ratios are summarized in Table [9l The systematic errors come mainly 
from selection criteria, event fitting, branching fraction uncertainties and re-weighting, while the 

3 To account for the asymmetry in 7r-A collisions, 
E791 used an additional offset parameter 
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-0.15 -0.125 -0.1 -0.075 -0.05 -0.025 0.025 0.05 

X F 

Fig. 6. Differential visible cross sections (—0.15 < if < 0.05) for D° and D + production: (a) da/dp^, with the fit 
of Eq.[T] and (b) da/dxp with the fit of Eq.[5]with a free parameter n and a fixed boundary parameter Xb = 0.062. 



Table 9. Ratios of cross sections. The first error is statistical and the second systematic. 



D+/D 0.41±0.06±0.04 

D* + /D° 0.44±0.11±0.05 

D+/(D+ + D°) 0.27±0.09±0.05 

D* + /D + 1.07±0.26±0.14 



luminosity error cancels. The value for the ratio a(D + )/a(D°) = 0.41 ± 0.06 ± 0.04 is the most 
accurate measurement of this ratio in pA reactions. It is in good agreement with the combined 
results from hadroproduction [lj as well as from e + e~ experiments [23]. The ratio also agrees 
with a simple prediction based on isospin symmetry and the measured ratio of vector to scalar 
meson production cross sections pQ. A comparison with results from other experimental studies 
is presented in Fig. [7J 

The ratio a(D* + )/a(D°) = 0.44 ± 0.11 ± 0.05 is also the most precise measurement of 
this ratio in pA reactions and is in good agreement with the results of NA27 and E769 [51 
[9]. The vector to scalar meson production ratio Py can be calculated in several ways (see [lj 
and references therein) if isospin invariance is assumed. From the ratio R± = a(D + )/a(D°) 
one obtains Py = (1 - i?i)/((l + Ri)Br(D* + -> D°tt+)) = 0.61 ± 0.09 ± 0.06. As a cross- 
check, we determine the same ratio from R2 = a(D* + )/a(D + ), and obtain the value Py = 
R 2 /(l + Br(L>* + -> D°tt + ) • R 2 ) = 0.62 ± 0.09 ± 0.05. The results are in good agreement with 
the world average value Py = 0.59 ± 0.01 pQ. 

Our result for the ratio a(D+)/(a(D°) +a(D + )) = 0.27±0.09±0.05 is the first measurement 
of this quantity in pA reactions. For comparison, the world average value of measurements in 
e + e~ collisions is 0.10±0.02 [24J, and 0.112+^0 in deep inelastic scattering at HERA [2]. 
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Fig. 7. Cross section ratio R = a(D + )/a(D°), comparison with previous experiments. The dotted line at R 
shows the prediction of the isospin model with Pv = 0.6 ,Tj; the solid line is a fit to the data points. 



0.42 



Table 10. Leading to non-leading particle asymmetries A in the visible range, —0.15 < ip < 0.05. The first error 
is statistical and the second systematic. 



D° 


D+ 


D*+ 


0.10±0.09±0.05 


0.14±0.11±0.06 


0.23±0.17±0.06 



5.4 Leading to non-leading particle asymmetries 



A leading particle is denned as one which has a light quark in common with the beam particle, 
in our case: anti-D mesons D , D~ and D*~ . The leading to non-leading particle asymmetry is 
defined as A = (olp — Cnon lp)/(clp + ^non lp)- The measured values of the asymmetry in the 
visible range, —0.15 < xp < 0.05, are given in Table [10] for the D°, D + and D* + mesons. The 
systematic uncertainty is dominated by the uncertainty of selection criteria and the kinematical 
re- weighting with respect to possible differences between leading and non-leading particles. Our 
measurements can be compared to the existing measurements of this asymmetry in different 
x F intervals: by E769 (0.06±0.13, 0.18±0.11 and 0.36±0.26 for D°, D+ and D*+ mesons with 
x F > 0.0 [9J), and E789 (0.02±0.06 for D° mesons with 0.0 < x F < 0.08 [8]). 



5.5 Atomic mass number dependence 



From the measured D meson production cross sections on three different target materials (Table 
[5]), the exponent a of Eq. 2] can be determined. The results of simultaneous maximum likelihood 
fits to the invariant mass distributions of individual material data samples for each D meson 
are summarized in Table [TTJ The systematic error has two contributions: the uncertainty in 
the luminosity per wire material (about 2.2%), and Monte Carlo statistics (less than 1%). The 
observed value, a = 0.99 ± 0.04 ± 0.03, is compatible with a linear dependence of the cross 
sections on atomic mass number (a = 1). Note that for the weighted average of a over all 
four samples, only those D* + events were considered for which the D° daughter particles were 
not reconstructed in the D° sample. Our result is in agreement with the result of E789 [5], 
a = 1.02 ±0.03 ±0.02. 
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Table 11. Atomic mass number dependence parameter a and its weighted average for the four D mesons. For 
the weighted average the value in parenthesis of the fifth row were used, corresponding to the subsample of D* + 
with D° daughter not common with the D° sample. 



Particle 


a 


D° 


0.969±0.057±0.026 


D+ 


1.051±0.082±0.028 


Df 


1.190±0.402±0.046 


D*+ 


0.832±0.138±0.022 




(0.847±0.185±0.022) 


Average 


0.994±0.044±0.025 



6 Summary 

With the HERA-B detector we have measured the total and single differential cross sections a, 
da jdp\ and da/dxF, the atomic mass number dependence of the cross sections, and the leading 
to non-leading particle asymmetries for the production of D°,D + ,Df and D* + mesons in pA 
collisions at the proton energy of 920 GeV. 

Extrapolating to the full phase space, the total cross sections per nucleon (in fib) are: 48.7 ± 
4.7 ± 6.6, 20.2 ± 2.2 ± 3.0, 18.5 ± 6.4 ± 4.1 and 21.6 ± 4.7 ± 3.6 for the D°, D + , D+ and D* + , 
respectively. In the range —0.15 < xp < 0.05 the measured cross sections are: 26.8 ± 2.6 ± 2.7, 
11.1 ± 1.2 ± 1.3, 10.2 ± 3.5 ± 2.0 and 11.9 ± 2.6 ± 1.7 for the D°, D+, D+ and D* + , respectively. 
The cross section per nucleon for cc production is ct(cc) = (49.1 ± 4.6 ± 7.4) ph. 

We have measured the cross section ratios a(D + )/a(D°) = 0.41±0.06±0.04 and a(D* + )/a(D°) = 
0.44 ± 0.11 ±0.05, as well as the vector to scalar meson production ratio, Py = 0.61 ±0.09 ±0.06. 
Our result for the ratio a(D+)/(a(D°) + a(D+)) = 0.27 ± 0.09 ± 0.05 is the first measurement 
of this quantity in pA reactions. 

From the measured atomic mass number dependence of the production cross section, the 
parameter a = 0.99 ± 0.04 ± 0.03 is extracted. This value is in agreement with the assumption 
of a linear dependence of cross sections, a=l. The measured leading to non-leading particle 
asymmetries in the xp range —0.15 < xp < 0.05 are consistent with existing measurements for 
different xf regions. 

The results of our studies are in good agreement with previous measurements of open charm 
production in pA interactions and provide, in the majority of cases, an improvement in accuracy. 
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